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APRESENTAÇÃO 


As ciências exatas, demonstram o espectro de conhecimento que 
desvenda os segredos do universo através da Matemática, da Física, da 
Química, da Computação e conexos. Um domínio onde a teoria e a prática se 
entrelaçam, no emaranhado que impulsiona o progresso tecnológico e industrial, 
por exemplo. E também tem conexão com a busca para teorizar o TODO, na 
medida do possível. 

De forma geral, as ciências exatas podem ser a base para o 
desenvolvimento de novas tecnologias que transformam a forma como vivemos, 
trabalhamos e nos relacionamos. A internet, por exemplo, é um produto da 
Física, da Computação e da Telecomunicação, enquanto os smartphones são 
resultados da convergência de diversas áreas, como Física, Química, Engenharia 
e Computação. 

Diante disso, o livro apresenta os conceitos teórico-práticos nos 
resultados obtidos pelos diversos autores e coautores no desenvolvimento de 
cada capítulo. Ademais, a Atena Editora oferece uma divulgação científica com 
qualidade e excelência, fundamental para garantir o destaque entre as melhores 
editoras do Brasil. 


Fabrício Moraes de Almeida 
Ciências Exatas, Computação e Engenharia. 
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RESUMO: O capítulo de livro apresenta 
fundamentos do estudo da previsão de 
demanda correlato a séries temporais 
de energia elétrica, desenvolvendo um 
processo de otimização para o atendimento 
no consumo de energia elétrica, com o 
objetivo de aprimorar a rotina de previsão. 
Além disso, é apresentada uma comparação 
de três modelos autorregressivos utilizados 
para a otimização mencionada, com o 
modelo de lógica fuzzy, dimensionado 
por meio da função de autocorrelação 
estendida. Os dados de consumo de 
energia elétrica apresentam uma estrutura 
de série temporal sazonal, suscitando um 
recorte histórico com as características de 
consumo de energia elétrica do Estado de 
Rondônia, e neste objeto foi implementada 


Capítulo 1 


uma metodologia de previsão por meio de modelos já sedimentados e em comparativo 
a um novo modelo, apresentado em Métodos de Identificação Fuzzy Para os Modelos 
Autorregressivos Sazonais Mediante a Função de Autocorrelação Estendida. Analisa-se, 
então, as performances dos modelos e aplicação para previsão de demanda de energia 
elétrica a curto prazo, 5 (cinco) dias úteis da semana, para fins de contratação de pacotes 
de energia elétrica junto as concessionárias distribuidoras de energia elétrica, obedecendo 
a legislação vigente sobre leilões e contratos de compra. No decorrer do trabalho, foram 
analisados os resultados de previsão de energia elétrica, pelos modelos apresentados, o 
modelo proposto em Sistema Fuzzy relacionado a Autocorrelação Estendida, sendo o mais 
satisfatório confirmado por erros de previsão em relação a demanda de energia elétrica para 
o Estado de Rondônia, e atendendo legislação sobre previsão e demanda de energia elétrica 
no Brasil. 

PALAVRAS-CHAVE: Previsão de demanda. Lógica Fuzzy. Função de Autocorrelação. 
Energia elétrica. 


ELECTRICITY DEMAND FORECAST USING FUZZY LOGIC AND EXTENDED 
AUTOCORRELATION FUNCTION 


ABSTRACT: The study of demand forecast correlated to time series of electric energy 
develops an optimization process for the supply of electricity, with the objective of improving 
the forecasting routine. In this study, it is presented a comparison of three autoregressive 
models used for the mentioned optimization, with the fuzzy logic model, dimensioned through 
the extended autocorrelation function. The electric energy consumption data presents a 
seasonal time series structure, provoking a historical cut with the electric power consumption 
characteristics of the State of Rondônia, and in this object was implemented a methodology 
of prediction by means of already sedimented models and in comparison to a new model, 
presented in Fuzzy Identification Methods for Autoregressive Models Using the Extended 
Autocorrelation Function. The performance of the models and application for short-term 
electricity demand forecasting, 5 (five) five business days of the week, was analyzed for the 
purpose of contracting electric energy packages with the electricity distribution concessionaires, 
in compliance with the legislation current agreement on auctions and purchase agreements. 
In the course of the work, we analyzed the results of electric energy prediction, by the models 
presented, the proposed model in Fuzzy System related to Extended Autocorrelation, being 
the most satisfactory one confirmed by errors of forecast in relation to the demand of electric 
power for the State of Rondônia, and complying with legislation on forecasting and demand 
of electric energy in Brazil. 

KEYWORDS: Demand prediction. Fuzzy Logic. Autocorrelation function. Electricity. 


INTRODUÇÃO 

O trabalho sobre análise de séries temporais com predição de sinal de energia 
elétrica tem sua aplicação na contratação de energia elétrica para distribuição pelas 
empresas concessionárias de energia elétrica. O processo de contratação é definido no 
Plano Mensal de Operação (PMO) de cada concessionária, e deve obedecer à legislação 
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no contexto da faixa de contratação de energia elétrica de 100% a 103% da demanda 
conforme ONS/2017 (BRASIL, 2017). A lógica fuzzy representa esta condição de forma 
matemática e através de relação conhecimento profissional e modelo, as várias séries ao 
longo do tempo podem ter um sinal representante em séries temporais são representadas 
por modelos matemáticos utilizados na previsão de demanda em curto, médio e longo prazo. 
Neste trabalho os modelos apresentados para efeito comparativo são os Autorregressivos 
de Média Móvel (ARMA), o Autorregressivo Integrado de Média Móvel(ARIMA), o Sazonal 
Autorregressivo Integrado de Média Móvel(SARIMA) e, por fim, o modelo estudado, Fuzzy 
Sazonal Autorregressivo Integrado de Média Móvel (FSARIMA) definido por Função de 
Auto Correlação Estendida (FACE). 


A LÓGICA FUZZY 


A Lógica Difusa ou Lógica Fuzzy é uma representação de um processo eurístico, na 
qual o conhecimento do especialista é utilizado para implementar melhorias e qualidade a 
um produto. A teoria de Conjuntos Fuzzy é definida como parte da fundamentação deste 
trabalho. A Figura 1 apresenta o sistema fuzzy de Mandani , definido como um processo 
no qual um sinal real ou crisp de entrada pertencente a um universo de discurso U é 
fuzzyficado, ocorrendo quando o número crisp é classificado em uma regra fuzzy através 
de um valor da função de pertinência. 

Em seguida, uma máquina de inferência faz a análise das regras válidas no banco 
de regras fuzzy, gerenciando-as conforme o processo em andamento, gerando um sinal 
de saída fuzzy, o qual será defuzzyficado, apresentando, ao final, um valor Crisp na saída, 
conforme atesta (WANG, 1999). A Lógica Fuzzy aplicada a séries temporais e utilizadas 
como processo de predição através de algorítimos utilizando a função de autocorrelação 
estendida (Carvalho Jr.; Costa Jr., 2017). A definição de conjuntos fuzzy é baseado em um 
dimensionamento do Universo de Discurso, definido por U, e sendo este um conjunto não 
vazio, um subconjunto de U, que representa as regras a serem utilizadas é definido como 
A, uma função de pertinência //tA, que representa um calor no intervalo de U > [0, 1], como 
demonstrado na Equação (1), onde x é o valor da variável a ser utilizado, e a representação 
matemática do conjunto A = ([x, yA(x; xe U). 
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Base de Regras Fzzy 


Fuzzyficador Defuzzyficador 


xeUu vev 
Conjuntos Maquina de Inferência Conjuntos 
fuzzy em U Fuzzy fuzzy em V 


Figura 1 — configuração básica de um sistema Fuzzy. 
Fonte: Adaptado de (WANG, 1999) 


0<r<I,sexeA; 
ala) = (1) 
nata 0,sez É A. 


O processo de fuzzyficação acontece quando uma variável de entrada x+ Crisp sendo 
classificado em um conjunto fuzzy AO, neste caso existe um valor de pertinência para xx. 
Caso existam sinais de ruído de fundo ou erros, no processo de fuzzyficação estes sinais 
os quais são indesejáveis, devem estes ser suprimidos de forma eficiente no processo. A 
simplificação dos cálculos é um passo importante para que o esforço computacional seja o 
mais brando possível ao longo do processo fuzzy. 

O fuzzyficador tipo triângulo define a relação de entrada do número Crisp x* mapeado 
no universo de discurso U e sendo definido em uma regra AO também relacionada a U, 
definindo uma função de pertinência conforme a Equação (2), onde bi é definido como o 
parâmetro positivo no processo e “x” é o operador t — norma, sendo utilizado como produto 
algébrico ou mínimo (Carvalho Jr.; Costa Jr., 2017), sendo definido o número fuzzificado. 


(1-8) ss (1- 268), selta = 2] <bi=1;2+,n. 


0, caso contrário 


fiadt) = (2) 


A descrição suppA é definida como o suporte do subconjunto A, conforme definido 
pela Equação (3) sendo um conjunto clássico, e sendo o processo de definição A como 
um subconjunto de U, e é definido por uma função de pertinência /A(x) no domínio de U e 
como contradomínio do intervalo [0, 1], conforme (WANG, 1999). 
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sppA=(xe U:u (x) >0) (3) 


Um conjunto fuzzy implementado sob uma restrição ou um limite estabelecido ao 
domínio, e considerando um valor a que é definido como um nível a ser aplicado as regras 
A do conjunto fuzzy, com a € [0;1], então o conjunto clássico A pode ser representado por 
[Aja=(xe U, yA > a). Sendo a> 0 e definindo [AJO como fecho do suporte 4, ou seja [AJO = 
suppaA, (BITAR et al., 2009). O a nível, indica o mínimo valor da função de pertinência uy A(x), 
e estes valores irão compor a saída do conjunto resultante fuzzy, conforme (CARVALHO 
JR., 2016). 

Quando as condições de atendimento, [AJO = 0 para0< a <1,[Alao qual é umintervalo 
fechado em R para O < a<1 eo subconjunto fuzzy A ter um suporte limitado, o valor da 
função de pertinência uy A(x) será definido como um número fuzzy. Uns dos procedimentos 
mais importantes para a lógica fuzzy são a determinação da função de pertinência, as 
escolhas dos modelos que representação as funções a serem representadas, conforme o 
especialista, uma função triangular, trapezoidal, gaussiana etc, conforme (ZADEH, 1965). 
Para um modelo triangular isósceles como modelo para regra 4, , definindo uma variável 
aleatória uma variável linguística Zi sendo i= 1,2,...,n,onde (n< T) e yA(x) um número 
fuzzy que representa o conjunto fuzzy triangular representado ela variável linguística Zi, 
conforme a Equação (4). 


0, se Zi< aj; 
Za) 
Z) = Rs se a;<2Z<b; 
gal à) = cj-Zi 
a, se b;<Z<c; 


0, se Zi; > C; 


Os parâmetros em sequência, aj, bj e cj, representam valores escalares indicando a 
posição dos vértices do triângulo isósceles e do ponto médio (valor de pertinência máxima). 
conforme a relação de Zi com as regras estabelecidas têm em seu resultado o valor da 
função f (Zi, aj; bj; c)) definidas nas Equações (4) e (5), na qual os vértices do triângulo 
apresentam aj< bj< cj sendo bjo pico do triângulo conforme descrito em (LAGO NETO, 
2011). 


f(Zua;byrc;) = mar [mim (pes; E) 0] (5) 
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As regras fuzzy que estão definidas como triangulares têm a anotação A = (aj; bj; 
cj), para a implementação do a-nível, (BITAR et al., 2009) apresentam as regras conforme 
a Equação (4.6),na qual a pode ser definido como demonstrado na Equação (7), a forma 
representativa do intervalo probabilístico ou de a nível na Equação (8). 


[4º =[a(b=a)+a;a(b,-c)+c] (6) 

log! a; Ee a; ev o = (7) 
b;-a; cb; 

[4º =<[aj"; [c])=[ a(b,- a) + a; a(c- b)+ cj (8) 


Além disso, a Lógica Fuzzy apresenta em sua álgebra conceitos como s-norma, que 
define as operações para união, e t-norma para operações de interseção, para a extensão 
do conjunto o operador complementar, conforme (WANG, 1999). A função de pertinência 
uA(x) e seu complemento é a função de pertinência u A(x) = 1 — /A(x), para a união de dois 
conjuntos do tipo 4 e B conforme a Figura 2 obtêm-se a seguinte Equação (9) e para a 
intersecção a Figura 3 esta relacionada a Equação (10). 


AU B 


B 


Figura 2 — Representação de União no sistema Fuzzy. 
Fonte: Adaptado de (WANG, 1999). 


Hu) = max [4 (9; u (X)] (9) 


Vans) = Min [pu (09; Gv0] (10) 
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Figura 3 - Representação de Intersecção no sistema Fuzzy. 
Fonte: Adaptado de (WANG, 1999) 


Para cada regra Aou B, devem estar contidos valores que representam um conjunto, 
exemplo uma série temporal, de um universo de discurso U, relacionados a função de 
pertinência de cada regras, dos valores de conjunto X pertencentes a U, os valores de 
entrada de parte real é definida como número Crisp, ou conjunto Crisp, conforme na 
literatura clássica. Os conjuntos que representam os valores de entrada são definidas em 
intervalos de números Crisp. A variável linguística representa estas faixas de valores, que 
são rotulados para cada conjunto fuzzy, definidos pelos especialistas, e o valor da função 
de pertinência referente a cada Regra. A variável linguística tem por definição o (X, T, U, 
e M) na qual X representa um termo que caracteriza um fenômeno, conceito, etc. .., To 
quanto é intenso ou pequeno com relação aos valores que X pode assumir, U descreve o 
universo que se enquadra o conceito de X, com valores Crisp de forma quantitativa tendo 
um valor mínimo e outro máximo, e M representa a semântica para a categoria de T dando 
valores linguísticos, como por exemplo “muito rápido”, “tensão baixa”, “alta potência”, 
relacionando X aos valores que assume em U, (WANG, 1999). 

As regras A, B, etc, que definem a Lógica Fuzzy, são compostas por termos como 
“SE”parte condicionante a uma faixa em referência, como “muito baixo”, “ângulo médio”; o 
termo “ENTÃO”"representa o resultado do acontecimento definido, denominados de regras 


Bar 


fuzzy “se - então”, conforme a Equação (11). 
SE < proposição fuzzy > ENTÃO < proposição fuzzy>. (11) 


O processo de defuzzyficação é o mapeamento de um conjunto fuzzy BO sobre o 
universo de discurso de saída Vc R sendo a saída da máquina de inferência conforme a 
Figura 4.1, sendo que o valor de saída y* é um número Crisp, e ys E V como apresenta 
(WANG, 1999), sendo o processo utilizado como defuzzyficador o Centro Médio ou 
centroide, conforme a Equação (2.12). O Defuzzyficador Centro Gravidade possui os 
critérios de plausibilidade o ponto y* pode representar o conjunto fuzzy BO do ponto de 
vista intuitivo, simplicidade computacional, por mais que seja um processo contínuo, a 
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representação é discreta em computação, não havendo problemas em continuidade com 
as variações pequenas no conjunto fuzzy B que não reflete em mudanças em y+, abrange 
todos as regras fuzzy B, com chances baixas de erro no resultado de saída do sistema e 
a implementação computacional é facilitada devido à utilização de regras baseadas em 
triângulos como regras fuzzy; e muito comum sua aplicação em sistemas de controle e 
como contribuição neste trabalho a séries temporais de energia elétrica (CARVALHO JR., 
2016). Segundo (WANG, 1999), sendo um processo computacional intuitivo e com uma 
entrada y,/BO(y) como uma relação do valor máximo em //BO(y) sendo o valor da função de 
pertinência para a saída e como resultado o número Crisp yx. 


is / ynp'(y)dy (19) 
v ue(y)dy 


A implementação do a como corte para as regras fuzzy BO é devido o cálculo na 
qual valores muito pequeno de ,/BO possam ser limitados, posto que os valores de ye V 
apresentam irrelevância para o cálculo de y+, sendo definida como uma saída defuzzyficada 
de centro de gravidade indexada por a conforme a Equação (2.13). Sendo que o universo 
de discurso indexado por a é dado por Va=(y e VIyBO(y) > a). 


e Jos yt'(y)dy (13) 
Jaz HB(y)dy 


Para todo sistema implementado há um processo, o ponto de segurança é 
fundamental, para este modelo de implementação há um processo, no qual o ponto de 
segurança é fundamental. Para este modelo de implementação, devemos estudar o 
intervalo de confiança fuzzy, Segundo (CARVALHO JR., 2016) as variáveis aleatórias 
independentes apresentam suas variações mediante um intervalo aleatório [L(X), U(X)], 
onde estas variáveis são identicamente distribuídas, sob intervalos onde X1, ---, Xn, e 
representado pelo vetor X= X7,---, Xn, L(X) > U(X) e deve satisfazer a igualdade conforme 
a Equação (14), e consonante a 6 como intervalo de confiança, os valores estimados fuzzy 
são identicamente distribuídos a série real e na Fígura 4. 
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Figura 4 — Representação de uma saída no Sistema Fuzzy através do método do centro de gravidade, 
aplicado ao intervalo de resultante na fuzzyficação das séries temporais de Energia Elétrica do Estado 
de Rondônia. 


Fonte: Adaptado de (WANG, 1999) 


PAUX) >0>U(X)=1=a (14) 


Baseado no inervado de confiança do sistema, em que Za representa o valor do grau 
de confiança onde o intervalo será indexado, o nível de significância a delimita gerando um 
nível cortando os valores das funções de pertinência e os valores abaixo deste nível não 
são computados, sendo o intervalo [L(x), U(x)], o intervalo de confiança para a variável 0 = 
1-a, onde os valores observados é definido por x=(xi), pertencente a variável aleatória X' 
onde i=1,2,:-:,n,e são vaiáveis aleatórias independentes e identicamente distribuídas e 
com fuzzy 67. A definição dos valores observados [" xjlLhe[” x]Uh,no qual he [0;1], por 
definição são sendo valores observados onde hX"iL h= X'ihLehXiUh=XihUi=1,2, 
***, n, e são variáveis independentes e identicamente distribuídas conforme apresentado 
em (CARVALHO JR., 2016). 

As Equações (15) e (16) representam o intervalos de confiança aproximados a (I.C) 
do modelo fuzzy, sendo os elementos shL = shu, que representam o desvio padrão para 
valores dos vértices do triângulo, e as termos xL he xUh, são definidos como a média 
entre os valores da base inferior e superior do triângulo 4, que define as regras do sistema 
fuzzy, com relação ao nível de significância a = 0, 05 e Za = 1, 96 para um intervalo de 
confiança de 95%, valores implementados em (CARVALHO JR., 2016). 


yn 
Le (ue) = [= (Cp + 289] = 1 (8) (6) (9 


Explorando o universo das ciências exatas: Teoria e aplicações Capítulo 1 


OS MODELOS APLICADOS A SÉRIES TEMPORAIS 


A série sazonal consiste em um conjunto de dados que em um período de tempo 
com início, meio e fim, suas características, repetem-se continuamente. A sazonalidade 
de um evento representa que este costuma ocorrer sempre num momento temporal 
específico, a aplicação as séries de energia elétrica é um exemplo quanto ao consumo 
de energia, apresentando um padrão que se repete de forma semelhante em um período 
de tempo. Quando há sazonalidade determinística, t é uma função determinística, e N é 
um processo estacionário, na Equação (4.44) obtêm-se o modelo sazonal determinístico 
com a aplicação da diferença sazonal. Quando há lags sazonais de múltiplos períodos, 
considera-se que exista uma sazonalidade estocástica, sendo este modelo dimensionado 
com os parâmetros p, d, qe P D, Qconforme a Equação (4.45), com o número de períodos 
sazonais s = 12, SARIMA (p, d, q) (P D, Q)12. 

Os dados observados ou minerados que definem um fenômeno físico e que produz 
um comportamento sazonal determinístico com período s=12, pode ser representado pela 
Equação (16), onde ut é uma função periódica que satisfaz o comportamento de yt= yt-12, 
onde outra relação é dada na Equação (17). 


Z=H+N, (16) 
(1-894,=0 (17) 


Onde Nt, sendo um processo estacionário e modelado por ARMA(p, q), apresenta- 
se como uma componente do sinal original, definindo a igualdade na Equação (18), sendo 
at uma função de ruído branco de média zero; porém a Equação (4.39) apresenta o ruido 
brando com sequência de senos e cossenos, onde constantes como uy, aj, Bj, onde j=1,2,:- 

, 6, não são definidas ou simplesmente não se tem conhecimento de seus valores, como 
solução é aplicado uma diferença sazonal conforme a Equação (20), e aplicando-se a 
Equação (17), conforme (MORETTIN; TOLOI, 2006). 


HBN= (Ba, (18) 


oe) (12xjt) 


=p 1008 cmi + À; 8EN n (19) 
J=1 E 
(1-B3)Z=(1-B3)y+(1-B3N, (20) 
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Em sequência de relacionamento aplicando-se (18) em (20), obtêm-se a Equação 
(22),e por adequação a Equação (23), onde : 


Ww,=(1- BZ, (21) 
(1-B9Z=(1-BNt (22) 
Bt - B%z, = 9 (Bt - Ba, (23) 
HB)W,= 0 (B1- B?) a, (24) 


A representação da série de energia elétrica de números convencionais reais(crisp) 
para um referente fuzzy e com saída novamente crisp, é a implementação do modelo de 
série temporal FSARIMA, (WANG, 1999). 

As séries fuzzy são implementadas neste processo analítico, com as entradas de 
números Crisp, o qual está contido no universo de discurso U. Seus valores são definidos sob 
regras, fuzzyficação; nestes casos representadas por triângulos, e saída, defuzzyficação, 
pelo método centroide, para uma resposta crisp, (ZADEH, 1965). As variáveis linguísticas 
que definem as regras geram, na função de pertinência PAi, valores definidos no intervalo 
[0,1], (WANG, 1999). 

De acordo com (SADAEI et al., 2014), o desenvolvimento da série temporal F (t, 
1) = F(t1, t-2) é definida como invariante em t, caso contrário a entrada F (?) será 
uma série de tempo variante distorcido, sendo F (t— 1) definida por Ai e F (?) por Aj sendo 
definido uma Relação Lógica Fuzzy (RLF), onde Ai relaciona Aj para à, j=1,2,3,... 
+ p, sendo p um subintervalo, definindo Ai o lado direito e 4j o lado esquerdo da ALF, 
as relações entre regras podem ser agrupamento em um grupo Difuso Lógico Ordenado 
(FLG), colocando todos juntos no RHS da FLG, e de acordo com (EFENDI; ISMAIL; DERIS, 
2015), podem ser escrito como Ai > Aj, Ai> Ak,..., Ai> Ap;i jk p,...,nmnonde n< 
T, para todo T igual a trajetória temporal de acordo com (CARVALHO JR., 2016). A serie 
temporal difusa F (7), que apresentam sazonalidade (FSARIMA) são des critos como uma 
relação de regras como deslocamento m, sendo F(t + m) > F(d), ou, F(t) = Ít + m), sendo t 
um período sazonal, (TSENG; TZENG, 2002). 

O valor que define o máximo valor da função de pertinência ,/Ai(x), relacionado a 
regra Ai, é um subconjunto do universo de discurso U, sendo definido por Supp(A) =(X 
L: 4Ai(x) > 0), sendo L um subconjunto difuso e X um conjunto compacto convexo, (BITAR 
et al., 2009). 

O sistema de referência a — cut é um conjunto variante no intervalo [0, 1] onde a 
representa o nível de referência para as regras [A] = (x € L : /A(x) > a), sendo a>0, e [AJO 
o encerramento do valor Supp(A), (WANG, 1999). 
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A implementação da Função de Autocorrelação Estendida Fuzzy (AECF) com a 
como conjunto de corte, é aplicado a função sazonal Wi,p,g a-cut, sendo i o grau do 
integrador, p do regressivo e q da média móvel, (CARVALHO JR., 2016). A ativação dos 
conjuntos contidos no universo de discurso L, serão relacionados a cada intervalo fuzzy, 
sendo os valores desta função os elementos da série temporal fuzzy com a-cut gerados 
pela a série original Zi(t), onde 1 <i<me 1 <i<T, sendo me T valores finais de cada regra 


fuzzy na Equação (4.61), onde p= 0,1,2,...,9=0,1,2,...,h=0,1,2,...,1,b=0,1, 
2,..., para pap-cuí(a, b) indica o operador com corte em especificado em a, (CARVALHO 
JR., 2016). 
ai E 
Wiza (h,0,b) = Zy"(a,b)—... 
ga, Zee" (a, b)-... (25) 


g-<i(a,b)Zee (a, b) —... 


Pp 


A FUNÇÃO DE AUTOCORRELAÇÃO ESTENDIDA. 


A função de autocorrelação estendida (face) foi proposta por (TSAY; TIAO, 1984) 
para definir de forma única os valores de p e q para dimensionamento do modelo ARMA(p,9) 
estacionário ou não. Conforme (MORETTIN; TOLOI, 2006), as funções de autocorrelação 
e autocorrelação parcial não são bons mecanismos para uma determinação precisa e com 
eficiência para definir os mesmos parâmetros p e q; em (CARVALHO JR., 2016) ela é 
utilizada para definir p e q em um sistema FSARIMA (p, d, q)(P, D, Q)12. Para a Equação 
(26) esta define como modelo ARMA estacionário para ser dimensionado pela face. 


Pj=dApja+...+ Okpj-k + 


+a(j)—...—Oa(j—1) (26) 


Como o índice j do processo pj é definido como j > 2, é definida uma matriz que 
determina os coeficientes dj (1), dj (2), dj (3) como no exemplo para 3 coeficientes da 
Equação (27). 


Pj Pj-1 Pj-2 od Pj+1 
(1) 
Pri Pj Pj- ds = |Pj+2 (27) 
Pj+2 Pj+1 Pj 69 Pj+3 


Define-se que para em qualquer posição de / na Equação (28), modelo sazonal, 
onde este deve ser maior ou igual ao regresso que possa zerar seu índice na Equação 
matricial conforme em a Equação (4.65),(CRYER; CHAN, 2008). 
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We =(1-0PB-..— BZ, (28) 


A Equação é redefinida para a Equação (29), 


Gióã) = Pl) (29) 


Sendo definido Gj, (k), P(jk) na forma de matriz e os coeficientes q&( (| k)) devem ser 
transpostos conforme a Equação (29), o número de elementos pk,j, na Equação (30), é 
definida como a função de autocorrelação estendida. 


do, * 
) Cen = 


et *(3) 
da) Crea) 28) 


Pkj = (30) 


Primeira análise, em 12, k = p, o processo Wk,t ()), é um processo de média móvel 
MA (9). 


| 

— +... + » j= 

mao da (1+6+... +64) to (31) 
0 «. J>p 


Quando k > p, é um processo estacionário ARMA (p, q), demonstrados nas Equações 
(31) e (32), (CRYER; CHAN, 2008). 


c, j=q+k-p 
= 32 
Ps | « j>qt+k-p (32) 


Onde -1<c<1. 
Para identificação de p e q no resultado gráfico do FACE é definido como a primeira 
sequência de zeros horizontal sem quebra de sequência e a diagonal em sequência de 


zeros sem interrupção, quando identificado o primeiro zero define p0 e g0 definindo estes 
valores conforme Tabela 1, conforme definida por (CRYER; CHAN, 2008). 
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q 

p 012... w-lqÔw Qt qt2.. 
0 EXXa XxX + X 

l LX] X 3% X 

2 LLX NL % 

ho xxx x 0 0/0 
ptl jxxxc xx 0/0 
pt? |jxxxa x x x 0 
ptd |Ixxxa xo x xo x 


Tabela 1 — Função de Autocorrelação Estendida simplificada, para dimensionamento do preditor de 
uma série temporal Fuzzy de energia elétrica aplicada a demanda no Estado de Rondônia. 


CONSIDERAÇÕES FINAIS 


Portanto, para os dados de consumo de energia elétrica presentam uma estrutura 
de série temporal sazonal, suscitando um recorte histórico com as características de 
consumo de energia elétrica do Estado de Rondônia, e neste objeto foi implementada 
uma metodologia de previsão por meio de modelos já sedimentados e em comparativo 
a um novo modelo, apresentado em Métodos de Identificação Fuzzy Para os Modelos 
Autorregressivos Sazonais Mediante a Função de Autocorrelação Estendida. 

Dessa foma, é possível uma análise das performances dos modelos e aplicação 
para previsão de demanda de energia elétrica a curto prazo, 5 (cinco) dias úteis da 
semana, para fins de contratação de pacotes de energia elétrica junto as concessionárias 
distribuidoras de energia elétrica, obedecendo a legislação vigente sobre leilões e contratos 
de compra. No decorrer do trabalho, foram analisados os resultados de previsão de energia 
elétrica, pelos modelos apresentados, o modelo proposto em Sistema Fuzzy relacionado a 
Autocorrelação Estendida, sendo o mais satisfatório confirmado por erros de previsão em 
relação a demanda de energia elétrica para o Estado de Rondônia, além disso, atende a 
legislação sobre previsão e demanda de energia elétrica no Brasil, podendo ser aplicado 
em outros tipos de demandas. 
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RESUMO: O satélite ICESAT-2 possui um 
sistema altimétrico que permite estimativas 
de superfícies. Este sistema possui um 
LIDAR na banda de 532 nanômetros que 
permite estimar altura de superfícies como 
icebergs, vegetação e também altura de 
superfícies aquáticas. Enquanto fótons 
emitidos com baixa energia permitem 
estimativas da altura de uma superfície 
aquática, fótons de maior energia permitem 
estimativas da batimetria destas superfícies. 
Dados do ICESAT-2 foram aplicados à 
região da Lagoa da Conceição, localizada 
na Ilha de Santa Catarina. Os resultados 
mostraram que estes dados permitem 
estimar estas alturas de superfície e 
batimetria para os dias particulares onde o 


Capítulo 2 


16 


satélite teve sua órbita sobre a região de estudo. Em particular, permitiu avaliar as regiões de 
maior e menor profundidades da lagoa. Os resultados mostraram que esta ferramenta pode 
ser uma importante opção para monitoramento da lagoa e suas possíveis mudanças em 
função de eventos meteorológicos extremos e tendências. 

PALAVRAS-CHAVE: ICESAT-2; Batimetria; Lagoa da Conceição 


APPLICATIONS OF THE ICESAT-2 SATELLITE FOR ESTIMATES OF SURFACE 
DEPTH FOR LAGOA DA CONCEIÇÃO 


ABSTRACT: The ICESAT-2 satellite has an altimetric system that allows surface estimations. 
This system has a LIDAR in the 532 nanometer band that allows estimating the height of 
surfaces such as icebergs, vegetation and also the height of water surfaces. While photons 
emitted with low energy allow estimates of the height of a water surface, higher energy 
photons allow estimates of the bathymetry of these surfaces. In this study, the ICESAT-2 
data were applied to the Lagoa da Conceição region, located on Santa Catarina Island. The 
results showed that these data allow estimating these surface heights and bathymetry for the 
particular days where the satellite had its orbit over the study region. In particular, it allowed 
the evaluation of regions having deeper and lower levels of the lagoon. The results show that 
this tool can be an important option for monitoring the lagoon and its possible changes due to 
extreme meteorological events and trends. 

KEYWORDS: ICESAT-2; Bathymetry; Lagoa da Conceição. 


INTRODUÇÃO 


A Lagoa da Conceição é uma lagoa tropical que está localizada na Ilha de Santa 
Catarina a leste do continente sul americano. A lagoa possui uma área aproximada de 
24 km? e um volume estimado de 66,6 x 10º mº* (AEROCONSULT, 2000). A lagoa possui 
um canal de ligação com o Oceano Atlântico que permite trocas de águas entre os dois 
sistemas. 

Medidas de batimetria da Lagoa da Conceição têm sido feitas no passado (Rudorff et 
al., 2005) e sua profundidade média estimada é de aproximadamente 1,7m a 8,9m (Godoy 
& Franco, 2007). No entanto, devido ao aporte constante de sedimentos, sua profundidade 
pode ser alterada periodicamente (Barbosa & José, 1998). 

Aregião costeira do sul do Brasil é constantemente afetada por marés de tempestade 
e marés meteorológicas (Melo Filho, 2017). A maré de tempestade pode ser definida como 
a provocada por ventos perpendiculares à costa que empilham as águas na interface com o 
continente (storm surges). Por outro lado, marés meteorológicas são variações do nível da 
água do oceano sobre a plataforma continental causadas por fenômenos meteorológicos. 
Estes fenômenos podem induzir tensões na superfície associadas à pressão atmosférica e 
ventos ao nível do mar que provocam variações no nível da água oceânica. 

Estudos recentes mostram que a interação do nível da água da lagoa com a do 
Oceano Atlântico funciona como um filtro de alta frequência, onde a maré astronômica 
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é filtrada pelo canal de ligação (Godoy et al., 2008; Bier, 2018). Portanto, estes estudos 
mostram que a maré meteorológica tem importante impacto na variabilidade da superfície 
da lagoa. Como as projeções climáticas apontam para uma maior frequência e intensidade 
dos fenômenos meteorológicos extremos (IPCC, 2021), torna-se de grande importância 
entender os mecanismos de interação e impacto dos fenômenos atmosféricos. 

O novo satélite Ice, Cloud and land Elevation Satellite-2 (ICESAT-2) colocado em 
órbita no final do ano de 2018, tem sido usado em várias análises, inclusive de altura de 
superfície e batimetria de lagos (Jasinski et al., 2024; Parrish et al., 2019; Markus et al., 
2017). Este satélite possui um laser a bordo que permite avaliar a altura de superfícies, 
como relevo, prédios, vegetação, e superfícies aquáticas. O satélite estima estas alturas 
a partir do tempo de ida e retorno dos fótons emitidos por seu laser. Além de superfície de 
águas continentais como lagos e rios, os dados do sensor permitem também estimativas de 
nível da água e da batimetria de águas oceânicas costeiras (Xu et al., 2022). Portanto, este 
sensor pode ser usado como uma importante ferramenta para monitoramento de níveis de 
superfícies aquáticas, como é o caso da Lagoa da Conceição. 

Este estudo teve como objetivo aplicar estas novas técnicas de análise por 
sensoriamento remoto para estimar a profundidade da superfície da água e da batimetria 
da Lagoa da Conceição por uso de dados do satélite ICESAT-2. 


METODOLOGIA 


Área de estudo 


A área geográfica abrangida por esse estudo consiste na região sobre a Lagoa 
da Conceição, localizada na região leste da Ilha de Santa Catarina do município de 
Florianópolis (Figura 01). 
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Figura 01 - Região de estudo, Ilha de Santa Catarina, Florianópolis e Lagoa da Conceição. Fonte: 
NASA, worldview, satélite sentinel-2. 


Satélite IceSat-2 


O satélite ICESat-2 possui um lidar a bordo que envia fótons para a superfície da 
Terra e assim estima as alturas da superfície a partir do tempo de viagem desses fótons 
de volta ao seu receptor (Newmann et al., 2019). Em suas aplicações, o satélite estima 
altura de superfícies como de relevo, vegetação florestal, gelo, ondas oceânicas, nuvens, 
camadas de aerossóis e superfícies aquáticas (Neuenschwander et al., 2020; Ramos-da- 
Silva et al., 2017). Além de superfície de águas continentais como lagos e rios, os dados do 
sensor permitem também estimativas de nível da água e da batimetria de águas oceânicas 
costeiras (Xu et al., 2022). 

O satélite ICESAT-2 possui uma órbita polar a 500 km de altura. Os dados de altura 
da superfície e batimetria são obtidos através do sensor ATLAS (Advanced Topographic 
Laser Altimeter System) a bordo do satélite. Este sistema de altímetro a laser usa um 
lidar na banda de 532 nm para medir o tempo de voo dos fótons e, portanto, a altura da 
superfície amostrada por seus seis feixes de laser (Neumann et al., 2019). Seus feixes de 
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fótons emitidos estão distantes 30 m dos três pares, que por sua vez estão espaçados de 3 
km (Fig. 2). O laser envia cerca de 20 trilhões de fótons por pulso e apenas cerca de uma 
dúzia retorna ao seu telescópio. Estimativas da altura podem ser obtidas comparando o 
dado com o elipsóide WGS84 (Neuenschwander et al., 2020). 


flight direction 


Figura 02 - Características dos feixes de fótons do satélite IceSat-2. Três pares de feixes são emitidos 
que possuem distância aproximada de 3 km. Fonte: NASA. 


Dados ICESAT-2 


Dados do satélite ICESAT-2 foram usados para estimar a profundidade do nível 
da água e da batimetria da Lagoa da Conceição. As passagens deste satélite de órbita 
polar pela região da lagoa ocorrem aproximadamente a cada 16 dias permitindo assim 
estimativas periódicas das alturas de superfície e batimetria. Os dados disponíveis foram 
usados para avaliar a variabilidade destas alturas. Estes dados de satélite são de livre 
acesso e podem ser obtidos do portal da NASA ou em https://nsidc.org/data/. 

O portal da NASA hitps://openaltimetry.earthdatacloud.nasa.gov/data/icesat2/ foi 
usado para acessar os dados e verificar a possível passagem sobre a região da Lagoa da 
Conceição. A Figura 03 apresenta um exemplo de estimativa da passagem do satélite com 
uso deste portal. 
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OPENALTIMETRY 


O arios 
O atio? 
O arios 
O ariio 


Figura 03 - Exemplo de análise do portal openaltimetry para download dos dados de altura da 
superfície com foco na Lagoa da Conceição. Fonte: NASA. 


A Tabela 01 apresenta as datas em que o satélite sobrevoou a região da superfície 
da Lagoa da Conceição. Estas datas correspondem apenas aos dias em que foi possível 
estimativas pelo satélite, pois em caso de cobertura de nuvens fica impossível aos fótons 
chegarem até a superfície. 


Data - passagem ICESAT-2 


17/01/2019 15/04/2020 13/07/2021 11/01/2022 


18/04/2019 14/10/2020 28/11/2021 


Tabela 01 - Datas de passagem do satélite ICESAT-2 sobre a região da Lagoa da Conceição. 


RESULTADOS E DISCUSSÃO 


A seguir são apresentados resultados de estimativas de altimetria da superfície da 
Lagoa da Conceição para vários casos em que houve a passagem do satélite. 


Altura de superfície e batimetria 


Inicialmente, após a pré-seleção dos casos através do portal openaltimetry foram 
feitas análises de algumas regiões sobre a Lagoa da Conceição. São apresentados vários 
casos baseados na Tabela 01. As figuras 4 a 11 apresentam os locais de passagem do 
satélite e estimativas amostradas nas regiões destacadas. Os resultados mostram que o 
satélite permite estimar as alturas superficiais da lagoa quando comparados com o elipsóide 
de referência. As várias regiões amostradas apresentam grande diversidade de alturas. Em 
várias ocasiões é possível também verificar a profundidade da batimetria da região. 
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Na passagem do dia 17/01/2019, que ocorreu mais próximo da interface a leste 
da lagoa, é possível notar alturas de profundidade da ordem de 1 metro (Fig. 4). Esta 
profundidade é mais rasa nestes locais. No entanto, nas passagens que ocorreram mais 
no centro da lagoa (e.g. dias 18/04/2019, 15/04/2020, 13/07/2021), as profundidades são 
da ordem de 2 metros (Figs. 5, 6, 8). Na passagem do dia 28/11/2021 é possível notar a 
transição das regiões de águas mais profundas para as mais rasas (Fig. 9). 

Na passagem do dia 11/01/2022 não foi possível avaliar a batimetria (Fig. 10). No 
entanto, uma amostragem para a região oceânica próxima mostra importantes aspectos 
da natureza ondulatória da superfície (Fig. 11). Este tipo de avaliação marítima pode ser 
importante, pois variações do nível oceânico próximo, pode afetar os fluxos de água entre 


a lagoa e o oceano. 


Altura (m) 


-27,566 -27,564 -27,562 -27,560 -27,558 


Latitude 


£ 


Figura 4 - Local de passagem do satélite ICESAT-2 e profundidades estimadas para o dia 17/01/2019. 
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Altura (m) 


4 


-27,555 


-27,570 -27,565 -27,560 


Latitude 


Figura 5 - Local de passagem do satélite ICESAT-2 e profundidades estimadas para o dia 18/04/2019. 


Altura (m) 


-27,5575 


-27,5700 -27,5675 -27,5650 -27,5625 -27,5600 
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Figura 6 - Local de passagem do satélite ICESAT-2 e profundidades estimadas para o dia 15/04/2020. 
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Figura 7 - Local de passagem do satélite ICESAT-2 e profundidades estimadas para o dia 14/10/2020. 


Altura (m) 
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Latitude 


Figura 8 - Local de passagem do satélite ICESAT-2 e profundidades estimadas para o dia 13/07/2021. 
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Altura 28/11/2021 


Altura (m) 


-27.570 


Latitude 


-27,560 


Figura 9 - Local de passagem do satélite ICESAT-2 e profundidades estimadas para o dia 28/11/2021. 
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Figura 10 - Local de passagem do satélite ICESAT-2 e profundidades estimadas para o dia 11/01/2022. 
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Return Signal Photons 
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Figura 11 - Local de passagem do satélite ICESAT-2 e profundidades estimadas para o dia 11/01/2022, 
sobre o oceano próximo. 


CONCLUSÕES E DISCUSSÃO 


Estimativas da profundidade da superfície e batimetria da região da Lagoa da 
Conceição foram estimadas com dados do altímetro do satélite ICESAT-2. Os resultados 
mostraram que estes dados permitem estimativas da altura da superfície da lagoa nos 
momentos da passagem do satélite. Regionalmente, os dados permitem avaliar a 
profundidade da superfície aquática e assim monitorar suas possíveis mudanças. No 
futuro, sugere-se que sejam feitas medidas in situ para assim ter uma melhor comparação 
com estes dados. Esta nova ferramenta pode tornar-se um importante instrumento de 
monitoramento desta lagoa, pois é de grande interesse da sociedade. 
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RESUMO: Com o avanço das emissões 
antropogênicas de gases de efeito estufa 
na atmosfera ocasionou um desbalanço 
na concentração de gases como o CO,. 
Logo, necessita-se de tecnologias que 
visem a redução destas emissões. Para tal, 
projetos de captura e armazenamento de 
CO, (CCS) em meio geológico são soluções 
fundamentais, proporcionando uma matriz 
energética sustentável. O Brasil é um 
potencial contribuinte neste tipo de projeto 
devido a suas vastas bacias sedimentares, 
sendo uma delas a Bacia do Parnaíba, 
localizada no Nordeste brasileiro, e que 
contém campos de gás como no Parque dos 
Gaviões. Assim sendo, o presente estudo 
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visa avaliar o potencial de aplicabilidade de 
implementação de tecnologias de CCS na 
Formação Poti através da comparação de 
parâmetros petrográficos como porosidade, 
permeabilidade e profundidade do 
reservatório. Adicionalmente, condições de 
seguridade também devem ser analisadas 
como o afeito selantes das rochas. Os 
resultados obtidos definem a Formação 
Poti como um potencial reservatório de 
armazenamento de CO, com os parâmetros 
comparados sendo  contundentemente 
satisfatórios. Concernente a rocha 
capeadora, os reservatórios da formação 
estudada são capeados por intrusões 
ígneas que impedem estruturalmente o 
vazamento de CO, e possuem potencial 
para aprisionamento por mineralização, 
devido a reatividade com o CO,. 
PALAVRAS-CHAVE:  CCS; Formação 
Poti; Armazenamento de COZ; Parque dos 
Gaviões 


ABSTRACT: The  advancement cf 
anthropogenic emissions of greenhouse 
gases into the atmosphere has caused an 
imbalance in the concentration of gases 
such as CO,. Therefore, technologies are 
needed to reduce these emissions. To this 
end, Carbon Capture and Storage (CCS) 


Capítulo 3 


28 


projects in geological environments are fundamental solutions, providing a sustainable energy 
matrix. Brazil is a potential contributor to this type of project due to its vast sedimentary basins, 
one of which is the Parnaíba Basin, located in the Brazilian Northeast, and which contains 
gas fields such as Parque dos Gaviões. Therefore, the present study aims to evaluate the 
potential applicability of implementing CCS technologies in the Poti Formation through the 
comparison of petrographic parameters such as porosity, permeability, and reservoir depth. 
Additionally, safety conditions must also be analyzed, such as the effect of rock sealants. 
The results obtained define the Poti Formation as a potential CO, storage reservoir with the 
compared parameters being overwhelmingly satisfactory. Concerning the capping rock, the 
reservoirs of the studied formation are capped by igneous intrusions that structurally prevent 
the leakage of CO, and have the potential for trapping by mineralization, due to reactivity with 
CO,. 

KEYWORDS: CCS; Poti Foramtion; CO2 storage; Parque dos Gaviões. 


INTRODUÇÃO 


Com o avanço da atividade humana, o aumento das emissões de CO, na 
atmosfera tornou-se uma problemática importante ao longo dos anos, devido ao impacto 
ambiental decorrentes de gases de efeito estufa (Moreno, 2013). Essas emissões provêm 
principalmente da energia de consumo e do fato que 85% da energia primária é fornecida 
por combustíveis fósseis (Orr, 2004), fazendo com que a redução das emissões de CO, 
constitui um grande desafio a ser solucionado. 

Com isso, uma das medidas para reduzir as consequências ocasionadas pelas as 
alterações climáticas e limitar os potenciais efeitos do aquecimento global, é redução das 
emissões de gases com efeito de estufa (GEE), com ênfase no dióxido de carbono (CO,) 
(Nema et al., 2012; Leung et al., 2014; Papadis e Tsatsonis, 2020). 

Dentre os possíveis meios para mitigar as emissões de GEE, o sequestro geológico 
de CO, se torna um candidato cada vez mais bem visto como uma alternativa viável para 
reduzir a libertação de gases com efeito de estufa na atmosfera. Com isso, a Captura e 
Armazenamento de Carbono (CCS, sigla em inglês) apresenta-se como uma ferramenta 
com potencial significativo para reduzir as quantidades de dióxido de carbono numa 
perspectiva da utilização prolongada de combustíveis fósseis nas próximas décadas (IPCC, 
2005). 

Para remover esse gás da atmosfera, realiza-se um projeto CCS, onde é realizada 
a captura de CO, de indústrias com elevadas emissões e a sua injeção em formações 
geológicas, tais como aquíferos e reservatórios de hidrocarbonetos esgotados. Porém, 
um dos principais obstáculos para o armazenamento de CO, a longo prazo e em grande 
volume em formações geológicas é a garantia de um armazenamento seguro e confiável 


desse gás (Moreno, 2013). 
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Nas últimas décadas, métodos para o armazenamento de CO, em larga escala 
foram identificados e amplamente utilizados. Dentre eles, os reservatórios de petróleo e 
gás esgotados são atualmente o tipo de reservatório mais utilizado no mundo, funcionando 
na estratégia de recuperação avançada de petróleo enquanto os reservatórios de carvão 
não minerável, que podem ter o CO, adsorvido na matriz de carvão e produzir metano. 
Além disso, outros tipos possíveis de reservatórios são os armazenamentos em cavernas 
de sal minadas e armazenamento oceânico e formações salinas profundas (Bachu, 2000). 
Outrossim, os aquíferos profundos também são utilizados para injeção de resíduos líquidos 
perigosos e não perigosos, pois contêm água conata fóssil e de alta salinidade que não 
é adequada para uso industrial e agrícola ou para consumo humano. Com isso, as altas 
pressões encontradas em aquíferos profundos indicam que estes podem suportar a injeção 
de CO, (Moreno, 2013). 

Entretanto, para que os projetos de CCS sejam instaurados para que se obtenha a 
notável redução destas emissões de CO, atmosférico, bilhões de toneladas ou mais devem 
ser anualmente estocadas com um aumento significativo de injeção deste gás em regiões 
que haja o aumento das emissões (Benson & Cole, 2008). 

Com isso, os projetos de CCS constituem em duas formas de aprisionamento, os 
quais são classificados na forma estrutural e na forma química (IPCC, 2005); com relação 
ao estrutural (estratigráfico ou residual), esta forma é responsável por grande parte da 
retenção do composto químico injetado de maneira similar a formação das reservas de 
hidrocarbonetos, enquanto que a forma química (solubilizada ou mineralizante) de certa 
forma é a mais desejável, devido a sua credibilidade em retenção por períodos geológicos 
sem efeitos danosos ao comportamento geomecânico do meio poroso (Rackley, 2010; Tao 
e Bryant, 2012). 

No entanto, estudos realizados no Brasil sobre a capacidade de armazenamento 
de CO, revelam o alto potencial de retenção deste gás nas bacias sedimentares do país, 
sobretudo na área de offshore (Ketzer et al., 2015; Goulart et al., 2020). A seleção das 
regiões brasileiras com elevados potenciais para aplicação da CCS é um aspecto decisivo 
em relação às emissões globais (Brasil, 2015). Com relação às tecnologias CCS no Brasil, 
o governo possui um retrospecto favorável na aplicação de tecnologias CCS pelo setor 
privado na implementação e uso (Camara, 2003); porém, existe uma desproporcionalidade 
de produção nas bacias do sudeste em comparação com o Nordeste, em parte devido à 
proximidade com o polo econômico mais importante do país e a falta de projetos de ação 
e em parte porque detém as maiores reservas nacionais de petróleo. Entretanto, pode-se 
utilizar campos em bacias menos visadas para lidar com emissões locais em outras regiões 
e implementar projetos-piloto nestes locais (Ciotta, 2021). 

Desta forma, o objetivo deste trabalho é realizar um estudo avaliando o potencial 
de aplicação de projetos CCS na formação sedimentar Poti, na bacia do Parnaíba, através 
da comparação entre parâmetros admissíveis para projetos de injeção de CO,, visando a 
estrutura porosa e a seguridade do reservatório. 
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GEOLOGIA LOCAL 


A Bacia do Parnaíba, antes denominada Bacia do Maranhão ou Maranhão-Piauí, 
fica localizada na Região Nordeste do Brasil e abrange uma área de mais de 600.000 
km (Cordani et al., 1984). Localizada nos estados do Maranhão e Piauí (Figura 1), os 
seus sedimentos cobrem embasamento maciços do período do pré-cambriano que seriam 
produtos da fissão do supercontinente Rodínia (Cordani et al 2009). A litologia da bacia é 
bem diversificada, varia de folhelhos, evaporitos, arenitos, diamictitos, etc. A deposição da 
Bacia se deu no final do Ordoviciano até o cretáceo, com variados ambientes de deposição, 
como glacio-fluvial, com possível diamictitos de origem glacial na Formação Ipu (Caputo, 
1984), desértico na Formação Samambaia, fluviolacustre da Formação Pasto bom, entre 
outros ambientes (Vaz et al, 2007). 
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Figura 1. Mapa de profundidade da base da Bacia Parnaíba. 


Fonte: Miranda et al. (2018) 
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Os reservatórios analisados neste estudo estão localizados no Parque dos Gaviões, 
uma região no estado do Maranhão que abriga 7 campos de gás natural. O Parque dos 
Gaviões possui um sistema petrolífero atípico. O Parnaíba passou por importantes eventos 
magmáticos ao longo de sua história, que foram responsáveis pela formação de diques 
e soleiras na bacia. Os dois principais eventos magmáticos toleíticos do Parnaíba são: 
Mosquito e o Sardinha (Trosdtorf et al, 2016). 

Os eventos de magmatismo descritos acima influenciaram diretamente na formação 
de hidrocarbonetos na Bacia do Parnaíba. Aplicado ao caso das rochas geradoras da 
Formação Pimenteiras, a espessura da Bacia não é o suficiente para que as rochas dessa 
formação sofram maturação e forme hidrocarbonetos, porém as intrusões magmáticas 
fizeram com que essas rochas sofressem maturação e começasse a formar hidrocarbonetos. 
Assim, como essas intrusões influenciaram a geração de hidrocarbonetos na bacia, as 
mesmas provocaram a formação de um sistema de diques e soleiras que servem como 


armadilhas para a cumulação de hidrocarbonetos (Figura 2) (Miranda et al, 2018). 


RES DiPaám BUPoEm  ilogêfm  ElCadeças Em MiPmerisras Fm Bltam Fm. 


Figura 2. Seção 2D aleatória com as principais formações que fazem parte da bacia, juntamente com o 
sistema de diques e soleiras (em vermelho). É possível observar essas intrusões dentro de formações 
de rochas geradoras, como os folhelhos da Formação Pimenteiras (azul) e na Formação Cabeças 
(laranja) e Poti (bege) consideradas como formações reservatórios, na qual serve como armadilha de 
hidrocarbonetos. 


Fonte: Miranda et al, 2018. 
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Os reservatórios estudados pertencem à Formação Poti, que é composta por 
arenitos cinzas esbranquiçados com lâmina de siltito e folhelho (Pimentel, 2020), a qual 
possui uma profundidade de aproximadamente 1100 m e sua espessura compreende 
aproximadamente 300 m. A formação é dividida em três zonas, utilizando a qualidade 
do reservatório como parâmetro: (1) zona superior apresenta uma boa qualidade de 
reservatório, com maior homogeneidade em sua litologia, (2) Zona intermediária apresenta 
uma média qualidade de reservatório e (3) Zona inferior que apresenta uma má qualidade 
de reservatório (Miranda et al, 2018). 

As propriedades dos reservatórios de arenitos da Formação Poti dos campos de gás 
encontrados na Bacia do Parnaíba, possuem porosidade média de 18%, podendo variar de 
5% a 27%, com uma permeabilidade de 240 mD, (Miranda et al, 2018). 

Entretanto, Pimentel (2020) fez interpretação litológica do poço 4-0GX-49-MA, a 
qual apontou que a porosidade varia de 9% à 21%. A ocorrência da unidade é no intervalo 
1163/1659m, e nesse local há uma intrusão de diabásio de aproximadamente 170m de 
espessura e a ocorrência de rochas metassedimentares em contato com o topo e base da 
intrusão. 

Algumas regiões da Formação Poti se encontram seladas pelos diques e soleiras 
resultantes das intrusões magmáticas. Um exemplo é um esquema da seção geológica 
feita por (Cunha et al., 2012) (Figura 3). 


Depth (m) 


Embasamento 


Figura 2. Esquema estratigráfico do Parque dos Gaviões. 
Fonte: Cunha et al., 2012. 
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O processo metodológico deste trabalho consiste em uma revisão acerca do 
arcabouço geológico a fim de detalhar os parâmetros petrofísicos da Formação Poti, 
pertencente a Bacia do Parnaíba, bem como a estruturação geológica das formações 
sedimentares e intrusões ígneas. A revisão é baseada em trabalhos mais recentes voltados 
para a área mencionada da bacia e, além disso, trabalhos tais que trazem vertentes 
geofísicas que mostram a capacidade de seguridade no armazenamento do CO, injetado. 

Para avaliar a aplicabilidade da implementação de um projeto de CCS na formação 
estudada serão comparados os dados obtidos para a Formação Poti com as condicionantes 
de seguridades comentadas por Ketzer (2012), Mineiro et al. (2022) e Silva et al. (2022), bem 
como com os parâmetros descritos como imprescindíveis para a aplicação determinados 


por Callas et al. (2022) e que podem ser observados na Tabela 1 a seguir. 


Descrição Critério 
Profundidade até o topo da formação 1000 - 2000 m 
Porosidade > 10% 
Permeabilidade >10mD 
Espessura do reservatório >10m 


Tabela 1 - Parâmetros mínimos para aplicabilidade de um projeto de CCS em reservatórios de meio 
poroso 


Fonte: Adaptado de Callas et al. (2022) 


Por fim, após a comparação com os referidos parâmetros e condicionantes, 
comentar-se-á sobre a viabilidade e potencialidade da implementação dos projetos de 
acordo com as evidências observadas. 


RESULTADOS E DISCUSSÃO 


A importância da profundidade do reservatório para a aplicação de projetos de 
armazenamento de CO, torna-se um fator criterioso, pois abaixo de 800 m de profundidade 
o CO, encontra-se em se estado supercrítico, o qual o gás comporta-se como fluido devido 
a densidade típica, apresentando, portanto, um menor volume (Bachu et al., 2000; IPCC, 
2005; Ketzer et al., 2016). 

Além disso, os reservatórios devem atender critérios de segurança e integridade 
para armazenamento eficiente de fluidos e com menores riscos ao meio ambiente. Portanto 
o reservatório deve apresentar capacidade economicamente viável, ter taxa de injetividade 
razoável e boa capacidade de retenção pela rocha capeadora para evitar escape de CO. 
Não obstante, espera-se um ambiente geoquímicamente e geomecanicamente estável 
(Ketzer et al., 2016; Mineiro et al., 2022). 


Explorando o universo das ciências exatas: Teoria e aplicações Capítulo 3 


34 


Logo, observando os parâmetros da Formação Poti, de forma geral, os valores 
se ajustam muito bem aos valores de referência aceitáveis para potencial implantação e 
desenvolvimento de um projeto de armazenamento de CO,, salvo a pouca representatividade 
de porosidade no conjunto da formação que atingem porosidade de 9% (Tabela 2), a 
qual limita-se pela inferioridade do valor ótimo estipulado por Callas et al. (2022) para 
reservatórios maduros. 


Tabela 2 - Parâmetros da Formação Poti para aplicabilidade de projeto CCS 


Descrição Formação Poti 
Profundidade até o topo da formação 1100 m 
Porosidade 9-21% 
Permeabilidade 240 mD 
Espessura do reservatório 300 m 


Concernente a permeabilidade e espessura do reservatório, por serem parâmetros de 
valores consideráveis, o potencial de aplicabilidade aumenta impactando até na capacidade 
de armazenamento, uma vez que proporciona uma taxa de injetividade maior. Porém, não 
implicando em maior seguridade, pois de acordo com Pawar et al. (2022), maiores plumas 
de CO, podem evidenciar maior potencial de vazamento de fluido armazenado, devido ao 
espalhamento e maior área de contato do CO, no reservatório. 

Adicionalmente, sobre a seguridade no armazenamento, de maneira positiva 
destaca-se que estratigraficamente a camada a ser submetida para injeção de CO, 
tem camadas selantes tanto acima do reservatório quanto a abaixo. A parte superior é 
limitada por intrusões ígneas na bacia sedimentar, o que pode ser benéfico, pois o CO, 
por reatividade pode mineralizar na região de contato. Ao passo que, na parte inferior uma 
camada fina de folhelho da Formação Longá, que sela a mobilidade descendente do CO, 
por convecção. 


CONCLUSÃO 


A avaliação preliminar do potencial de aplicabilidade para operação de um futuro 
projeto de armazenamento de CO, por meio de levantamentos geofísicos e petrofísica 
de certa forma é um balizador determinante para o gerenciamento e futuras tomadas de 
decisões em projetos deste tipo. 

Logo, a Formação Poti atende às exigências mínimas interpostas para a 
implementação de um projeto de CCS, uma vez que, como comentado anteriormente, 
comparativamente seus parâmetros se enquadram perfeitamente além de já possuir 
instalações físicas como poços, proximidade de fontes emissoras e excepcionalmente 
possuir boa capacidade capeadora superior e inferior demarcadas pela soleira e pelo 
folhelho da formação Longá. 
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Deve-se, portanto, avaliar futuramente por meio de simulações de fluxo de fluidos 
a Formação aqui estudada com o intuito de averiguar a capacidade de armazenamento, o 
potencial de retenção, os efeitos geoquímicos de trapeamento ocasionados pela capeadora 
ser uma rocha basáltica. 
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RESUMO: O primeiro Geopark do Continente Americano, foi criado em 2006 - Geopark 
Araripe, que abrange uma área de 3.441 km? em torno de algumas cidades do Ceará, como 
Crato, Juazeiro do Norte, Barbalha, Nova Olinda, Santana e Missão Velha. O Geossítio 
Floresta Petrificada do Cariri, distante de Missão Velha 6 km, cujo acesso é por uma rodovia, 
trata-se de uma área privada cercada e sua entrada é controlada por cadeado. Possui uma 
carência de infraestrutura de apoio, pois possui grande potencialidade para o turismo. Nesta 
perspectiva, o objetivo geral deste estudo é compreender as formações dos troncos fósseis 
silicificados da formação Missão Velha e relatar as vantagens para o turismo e as dificuldades 
de implantação no Geossítio Floresta Petrificada. Adotou-se como metodologia uma Revisão 
Bibliográfica, desenvolvida a partir da junção dos trabalhos já publicados. Utilizou-se como 
plataformas de buscas: SCIELO, Google acadêmico, EDU e Periódicos da CAPES. Na 
busca por estudos relevantes, foi adotado “AND” como operador booleanos, fazendo uso 
dos seguintes descritores: Troncos fósseis “AND” Coníferas, Formação Missão Velha “AND” 
Jurássico. Os resultados foram desenvolvidos com base em 33 trabalhos, após adoção dos 
critérios de inclusão e exclusão. Os estudos nos mostraram que este geossítio possui uma 
exuberância acentuada em virtude da presença dos troncos fósseis, bem como um amplo 
potencial didático, associado à sua relevância para Formação Missão Velha. Portanto, este 
trabalho foi de total relevância para produção de conhecimentos acerca desta temática, no 
entanto, apresentou como principais limitações a restrição quanto os acessos de artigos 
publicados de forma não pública. 

PALAVRAS-CHAVE: Troncos Fósseis, Sedimentação, Fossilização, Fragmentos. 


GETTING TO KNOW THE PETRIFIED FOREST OF MISSÃO VELHA: A 
LITERATURE REVIEW 

ABSTRACT: The first Geopark on the American Continent was created in 2006 - Geopark 
Araripe, which covers an area of 3,441 km? around some cities in Ceará, such as Crato, 
Juazeiro do Norte, Barbalha, Nova Olinda, Santana and Missão Velha. The Petrified Forest of 
Cariri Geosite, 6 km away from Missão Velha, which access is via a road, is a private fenced 
area and its entrance is controlled by a padlock. It has a lack of support infrastructure, as it 
has great potential for tourism. In this perspective, the general objective of this study is to 
understand the formations of the silicified fossil trunks of the Missão Velha formation and to 
report the advantages for tourism and the difficulties of implantation in the Petrified Forest 
Geosite. A Bibliographic Review was adopted as a methodology, developed from the junction 
of already published works. It was used as search platforms: SCIELO, Google academic, 
EDU and Periodicals of CAPES. In the search for relevant studies, “AND” was adopted as 
a Boolean operator, making use of the following descriptors: Fossil trunks “AND” Conifers, 
Missão Velha Formation “AND” Jurassic. The results were developed based on 33 studies, 
after adopting the inclusion and exclusion criteria. The studies showed us that this geosite has 
an accentuated exuberance due to the presence of fossil trunks, as well as a wide didactic 
potential, associated with its relevance for the Missão Velha Formation. Therefore, this work 
was of total relevance for the production of knowledge about this theme, however, it presented 
as main limitations the restriction regarding the access of articles published in a non-public 
way. 

KEYWORDS: Fossil trunks, Sedimentation, Fossilization, Fragments. 
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INTRODUÇÃO 


Os Geoparques possuem condutas indispensáveis para conservação da natureza 
e aspectos abióticos, pelo universo, totalizam 161 geoparques em 44 países, na qual 
constituem o programa global de Geoparques da Unesco (Unesco Global Geoparks, UGG), 
o que identificou um aumento de medidas sustentáveis do meio ambiente, com notável 
cautela as comunidades que abrangem a área (NASCIMENTO et al., 2020). 

Segundo Brilha (2009), o Geoparque é um território com delimitações para a 
unidade, e toda esta área está voltada para fomentar o turismo e aumentar a produtividade, 
gerando uma forma de renda para os habitantes locais, fundamentado na preservação 
das riquezas paleontológicas, agregadas com os bens natural e cultural que formam o 
patrimônio. Macedo e Pinheiro (2014) enfatizam que o Geopark é uma área de grande 
relevância geológica, comprovado pelo seu encanto. 

O Geopark Araripe, abrange parte territorial de 3.441 km? em torno de algumas 
cidades do estado do Ceará, incluindo o Crato, Juazeiro do Norte, Barbalha, Missão Velha, 
Nova Olinda e Santana do Cariri (MOCHIUTTI et al., 2012). A Bacia do Araripe, é uma Bacia 
no interior do Nordeste, que está localizada no Sul do estado do Ceará, sendo a maior 
e mais completa, e com extensões que se limitam pelo estado do Piauí e Pernambuco 
(FREITAS et al., 2008). Abrange uma área de 12.000 km? (SARAIVA et al., 2007). A Bacia 
do Araripe ocupa um espaço em torno de 9000 km2, compreendendo os estados do Ceará, 
Piauí e Pernambuco (CHAGAS, 2006). 

O Sitio Ribeirão na cidade de Brejo Santo - CE, abriga fragmentos de madeira fóssil 
de grandes dimensões, possuindo de 160 a 220 cm de comprimento por 90 a 110 em 
de diâmetro (FREITAS et al., 2008). Outros locais aonde houveram pesquisas das quais 
descreveram fragmentos com mais de 1 m de comprimento, foram os sítios Várzea Grande 
e Pau Branco, na cidade de Mauriti - CE. Nestas descrições estudadas, as toras fossilizadas 
exibem somente o xilema secundário preservado, não existindo o córtex nem a medula. 

Na cidade de Missão Velha temos o Geossítio Floresta Petrificada, situado em um 
local chamado de Grota funda. Neste local é possível encontrar a melhor exposição natural 
da Formação Missão Velha, em toda a Bacia do Araripe, podendo servir de base para 
compreensão de toda a sequência sedimentar constituída (FREITAS et al., 2008). Este 
Geossítio abriga grandes paredões rochosos, onde a estratigrafia do terreno pode ser bem 
observada, bem como a formação Missão Velha (FAMBRINI et al., 2011). 

Mesmo diante da infinidade de trabalhos que discorrem sobre a geologia da Bacia 
do Araripe, percebido, nas últimas décadas, a ausência de argumentos mais precisos que 
definam a formação Missão Velha (FREITAS et al., 2008). Compreendendo que a falta de 
uniformidade de opiniões quando citam a estratigrafia da formação Missão Velha, faz surgir 
várias propostas argumentando sobre a proposição estratigráfica e conteúdo litológico 
sobre o que seria a formação Missão Velha. Neste sentido, a falta de uma explicação mais 
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proveitosa, está no fato da ausência de fósseis que orientem nesta discussão (FAMBRINI 
etal., 2011). 

Freitas et al. (2008) enfatiza que há 45 anos já havia informações sobre estas 
madeiras fósseis silicificadas, no entanto, existe uma carência de respostas envolvendo a 
identificação de sua tafonomia e adequando em um contexto geológico, que intensifique a 
necessidade de um aprofundamento desta área fossilífera. 

Desde o século XX, já se tinha conhecimentos sobre as madeiras fossilizadas, 
porém essas informações foram insuficientes, apresentando sua visibilidade somente após 
a criação do Geopark Araripe, reconhecido em 2006 pela UNESCO (BARRETO JUNIOR 
et al., 2021). 

Apesar de sua relevância, falta investimentos e um olhar mais preciso, pelo número 
de pessoas que visitam e fazem suas pesquisas, e precisam de uma estrutura de apoio. 
Faz-se necessário a preservação desta área por parte dos órgãos competentes, buscando 
aumentar as fiscalizações, por ter uma ótima localização, próximo à rodovia aumenta 
radicalmente o contrabando destes itens de um valor inestimável para os pesquisadores 
(FREITAS, 2019). Em decorrência da exuberância visível de fósseis, criaram o Geossítio 
Floresta Petrificada do Cariri, reconhecida pela sua dimensão de madeiras fossilizadas e 
pela sua seção tipo, formação Missão Velha (FREITAS et al., 2008; SILVEIRA et al., 2012). 


MATERIAIS E MÉTODOS 


Tipo de estudo 


Esta pesquisa foi desenvolvida por meio de um estudo bibliográfico, com abordagem 
qualitativa e objetivo exploratório. De acordo com Gil (2008), o objetivo de uma pesquisa 
exploratória é estar por dentro de todo o assunto sempre buscando mais estudo sobre o 
conteúdo como também ser, mas explorado sempre buscando informações para aprimorar 
o conhecimento. 

Já a revisão bibliográfica, é desenvolvido por meio da junção de todo material 
diversificado, já elaborado por autores que através de suas pesquisas contribuem para 
elaboração de posteriores estudos, as fontes de dados podem ser podem ser: livros, artigos 
científicos, periódicos de revistas, teses, dissertações entre outros (GIL, 2002). 


Instrumento e procedimento 


Foram selecionados 33 trabalhos científicos, com dados dos anos de 1983 a 2022. 
Utilizou-se como plataforma de estudo e buscas o Google acadêmico, o SCIELO, periódicos 
da CAPES e EDU. Para tanto, utilizou-se dos seguintes descritores: Troncos fósseis “AND” 
Coníferas, Formação Missão Velha “AND” Jurássico. 
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Critérios de inclusão e exclusão 


Adotou-se como critérios de inclusão, arquivos disponíveis de forma gratuita e 
integral nas plataformas supracitadas, que compreendessem os anos de 1983 a 2022, em 
língua portuguesa e que atendessem aos descritores de pesquisas. Os critérios de exclusão 
se voltaram para arquivos disponíveis de forma parcial, pago, e que não abordassem os 
descritores da pesquisa, e que estivesse faixa temporal definida. 


Análise de dados 


Os dados são discutidos e expressos de forma dissertativa, onde são expostas as 
principais informações que somam ao desenvolvimento deste estudo. 


RESULTADOS E DISCUSSÃO 


A Bacia do Araripe foi construída através de ocorrências tectônicas, determinante 
da fragmentação do Gondwania, ligada a abertura do Oceano Atlântico Sul, e o riftiamento, 
gerado através do afastamento das placas continentais da América do Sul e da África, 
que motivaram para a ocorrência da sedimentação e progressão desta bacia. No entanto, 
é mencionado um desentendimento quando é comentado a cronologia das rochas 
sedimentares basais da Bacia no que se refere a Formação Mauriti ou Tacaratu, que foram 
alvos de muita polêmica, assim como a formação Missão Velha (KELLNER, 2002). 

A formação Missão Velha fecha e inicia o clímax de Rifti e tem uma discordância 
interna. Próximo a Abaiara temos o Moro do Cruzeiro, não possui um nível de paleosolo, 
possui uma discordância erosiva apontado pela presença de fragmentos de arenitos 
inferior dentro do arenito superior. De modo que arenitos grossos com estratificações 
cruzadas acanaladas que são a passagem da sequência 1 para a sequência 2, isso seria 
o significado tectônico da Formação Missão Velha, isso acontece por que ela marca a 
passagem do Rifti, para o nosso entendimento. A grande discussão envolve este ponto. 
Em suma, o que acontece é uma Formação, mas, ela tem uma tecno sequência interna 
(FAMBRINI et al., 2019). 

A sequência clímax de Rifti neocomiano, compreende o item superior a formação 
Missão Velha e a formação Abaiara. Esta sequência evidencia uma desconformidade 
intraformacional onde é possível diagnosticar na formação Missão Velha, conglomerados 
de seixos areníticos com estratificação cruzada e acanalada e tabular, sabendo que podem 
se situar dispostos justamente por contato erosivo (FAMBRINI et al., 2011). 

Estes conglomerados partem literalmente e verticalmente para arenitos 
conglomeráticos amalgamados grosseiros, com estratificação cruzada acanalada e tabular 
de média escala, da qual podemos comprovar diversos fragmentos de madeira silicificada 
do gênero Dadoxylum benderii da qual chegaram a uma altura de 1,60 m na parte mais 
grossa da Formação Missão Velha (FAMBRINI et al., 2019). 
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Percebe-se uma ausência de conformidade quando se fala sobre as proposições 
estratigráficas da Bacia do Araripe, o que torna a coluna estratigráfica desta Bacia algo 
controverso, um exemplo disso é a unidade Missão Velha que expõe inúmeras sugestões 
estratigráficas que se divergem, quanto a posição, espessura, abrangência e conteúdo 
litológico. 

Estes conceitos que foram surgindo de vários autores acerca da Formação Missão 
Velha, não chegavam a um consenso definido para a unidade, devido à ausência de fósseis, 
que contribuam para uma idade confiável (FAMBRINI et al., 2011). A Bacia do Araripe exibe 
uma sucessão estratigráfica muito complexa, compostas pelas Formações Cariri, Brejo 
Santo, Missão Velha, Abaiara, Barbalha, Santana e Exu (CHAGAS, 2006). 

Assine (2007) e Fambrini et al. (2017) chegaram a um consenso sobre a coluna 
estratigráfica da Bacia do Araripe quando fazem a seguinte descrição: fracionaram em 
5 seguimentos sendo, a sequência Paleozóica representada diante a Formação Cariri, 
onde a sua composição tem arenitos médios a grossos fluviais de idade neordoviciana e 
eossiluriana. Sequência início de Rifti, neojurássica formada pelos pelitos da Formação 
Brejo Santo e também os arenitos da parte inferior da Formação Missão Velha. Sequência 
clímax de Rifti composta pelos arenitos grosseiros da parte superior da Formação Missão 
Velha e acrescentando os pelitos e arenitos da Formação Abaiara Neocomiana. Sequência 
pós-Rifti fracionada em duas, pós-Rifti |, Aptiano-Albiano, possuindo a seguinte composição 
as Formações Barbalha, Crato, Ipubi e Romualdo, e pós-Rifti II, Albiana-Cenomaniana, que 
se entende que a sua descrição é de sedimentos aluviais da Formação Araripina. 

Atualmente, sabe-se que existe deficiências de informações sobre a Formação 
Missão Velha, e isso requer, mas estudos comprobatórios desta formação, ressaltando 
a relevância na análise de fácies, reconhecendo a importância de novas pesquisas e 
profissionais interessados na busca por explicação para as seções colunares. 

Sabendo que a Bacia necessita de estudos estratigráficos com ênfase na observação 
de fácies de cada uma dessas unidades (CHAGAS, 2006). Pode-se afirmar, que a Formação 
Missão Velha foi sinônimo de muita polêmica por motivos de uma indefinição da unidade 
(ARAI, 2006). 

Como já frisado, a Formação Missão velha está sobreposta em concordância 
com a Formação Brejo Santo, possuindo dimensão que não ultrapassa o seu limite de 
200m, havendo composição arenitos quarizosos ocasionalmente feldspático, as vezes 
conglomeráticos dispondo de uma copiosa abundancia de madeira silicificada designada de 
uma floresta de coníferas a Dadoxilon benderi (ASSINE, 2007). É uma unidade, composta 
por arenitos finos, a grossos até conglomeráticos e isso evidencia pela presença de troncos 
fossilizados (FAMBRINI et al., 2017; ASSINE, 2007). 

Braun (1966) fundamentou a funcionalidade do termo, Formação Aliança. Diante dos 
atributos litológicos e as toras que se fossilizaram pertencente a Formação Missão Velha 
que podem ser encontradas com frequência a unidade, foi conferida como a Formação 
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Sergi, Dom João (BRITO, 1987). Fambrini et al. (2017) explica que a Formação Missão 
Velha foi diagnosticada com uma discordância erosiva interna, diferenciada pela existência 
notada dos arenitos conglomeráticos e conglomerados, bem significativos e atuantes em 
clastos de arenitos médios a finos, com procedência das rochas justamente oriundas 
abaixo da superfície. 

De acordo com Fambrini et al. (2017), a Formação Missão Velha passou por uma 
sequência de processos que corresponde: sequência 1 - possui abundância em arenitos 
médios e médias a finos, já a sequência 2, foi diagnóstica com seus arenitos grossos a 
muito grosso e de uma forma pouco recorrente encontraram conglomeráticos subordinados 
com uma quantidade bastante considerável de madeiras fósseis silicificados. Neste sentido, 
a Formação Brejo Santo representa a primeira fase lacustre da Bacia fornecendo indícios 
de um ambiente lacustre raso com influência fluvial. 

Foi verificado no Geossítio Floresta Petrificada do Cariri, situado numa localidade 
Rural, no qual o solo é produto da pedogênese, bastante desfrutado para a agropecuária, 
o que denota, como o solo foi reagindo a fatores (NASCIMENTO et al., 2020). Há autores 
que mesmo comprovando seus aspectos significativos explicam que não há necessidade 
de cautela, quando se fala em preservação da seção, mesmo diante de sua relevância não 
precisa de uma preocupação com preservação desta unidade (FREITAS, 2019). 

No entanto, alguns autores descrevem que o principal motivo de criação do Geossítio 
Floresta Petrificada do Cariri no Geopark Araripe, no ano de 2006, foi a remoção sem 
autorização destes sedimentos, o que motiva a uma perda simbólica destes sedimentos, 
mesmo compreendendo que esta informação é centenária da descoberta destas madeiras 
silicificadas, ainda há uma evacuação mesmo sendo de uma maneira parco (BARRETO 
JUNIOR, 2015). 

A Floresta Petrificada é um dos Geossítios do Geopark Araripe, contudo possui 
uma especificidade geológica muito limitada, tendo o número de visitas de discentes 
bem reduzido. Percebe-se ainda a falta de infraestrutura e apoio de órgãos competentes 
(GABRIELLI, 2014). No entanto, por ter só rochas areníticas e a ocorrências dos troncos 
silicificados, geralmente as pessoas que vão visitar, são, mas objetivas, ao acesso às 
informações cientificas. Sendo assim, por conta de sua grande especificidade geológica, 
leva a um reduzido número de visitação por público não acadêmico (FREITAS, 2019). 
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Ni 


Figura 1 — Floresta Petrificada de Missão Velha — CE 
Fonte: Autores (2023) 


Este geossítio não possui uma exuberância diversificada como encontramos na 
Formação Santana, onde é prestigiada por uma incrível quantidade de fósseis de peixes, 
insetos e pterossauros (BARRETO JUNIOR, 2015). No entanto, podemos encontrar uma 
paisagem exuberante, que em períodos chuvosos o verde proporciona uma paisagem 
encantadora da Grota funda que é seu cânion, e possui também várias trilhas, bem os 
paredões rochosos (FREITAS, 2019). 

Entre seus pontos negativos encontramos uma deficiência nos serviços fundamentais, 
que ajudarão estudantes e pesquisadores em obterem um melhor aproveitamento em suas 
pesquisas de campo. É importante fiscalização e manutenção das trilhas, pois caso estejam 
mal sinalizadas, isso pode fazer com que os visitantes se enganem na localização, levando- 
os a se perderem nas trilhas. Reforça-se então a carência de serviços fundamentais como 
ausência de infraestrutura de apoio, falta de energia, de abastecimento de água, trilhas mal 
sinalizadas (FREITAS, 2019). 

As toras que se fossilizaram durante o período geológico, podem ser visualizadas 
no Geossítio Floresta Petrificada do Cariri, são somadas a indícios do qual no período 
jurássico ocorriam extensas florestas de coníferas, totalmente distinta dos ecossistemas 
encontrados recentemente (NASCIMENTO et al., 2020). Os troncos fósseis caracterizados 
como Dadoxylon há pontos que supostamente levem a fazer parte das famílias Taxodiáceas 
ou Pináceas (FREITAS et al., 2008). 
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Devido a uma similaridade com troncos de Queirolepiácea identificadas na 
formação Navajo, Estados Unidos (STOKES, 1991). Os xilemas secundários da formação 
Missão Velha e formação Navajo possuem esta semelhança, no entanto, não é possível 
diagnosticar em detalhes a constituição de suas flores presentes no leste Cariri Cearense, 
Bacia do Araripe em períodos remotos ou que seus fragmentos estejam em ótimo estado 
de preservação (FREITAS et al., 2008). 

Quando se fala sobre o gênero das coníferas, Freitas et al. (2008), Freitas (2019) e 
Assine (2007) enfatizam que os copiosos troncos silicificados da Formação Missão Velha, 
sugere que sejam de Coníferas do gênero Daxoxylon Benderi. No entanto, não chegando 
a resultados conclusivos sobre sua família se é Taxodiácea, Pinácea ou até mesmo 
Queirolepiácea, por conta da falta de informação da sua medula e xilema primário. 


CONCLUSÕES 


A região do Cariri dispõe de uma grande geodiversidade, na qual se destaque os 
geossítios do Geopark Araripe, possibilitando diversas maneiras de realizar o geoturismo. 
Propostas relacionadas a geoconservação, são de extrema relevância, bem como a 
compreensão de formas de proteção e conservação deste patrimônio tão relevante. 

Diante das evidências, foi possível observar que falta um consenso por parte dos 
autores, quando é mencionado as proposições estratigráficas da Bacia do Araripe, um 
exemplo, é a Formação Missão Velha que já foi motivo de muita discussão. Em virtude de 
não possuir fósseis que contribuam para uma idade confiável. Trata-se de uma unidade 
onde são escassos os estudos científicos. 

O Geossitio Floresta Petrificada do Cariri, possui uma grande especificidade 
geológica, onde podemos evidenciar, a paisagem do Canio da Grota Funda, suas trilhas 
e os paredões rochosos que exibem a Formação Missão Velha, além da exuberância dos 
troncos fósseis. 

É possível encontrar exposições dos fragmentos da Formação Missão Velha 
distribuídos pelos municípios de Abaiara, Missão Velha, Mauriti e Brejo Santo. Nas 
localidades do Morro do cruzeiro, Serra do Espia, Linhas de Alta Tensão, Cachoeira do 
Macedo, expostos pelas formas in situ, ex situ e exumados, de uma assembleia Alóctones 
do gênero Dadoxylon benderi, precisando ainda de mas estudos sobre sua tafonomia e 
gênero. 

Desta forma, é possível concluir que há a predominância de arenitos amalgamados, 
verificando que a posição da qual estão os troncos fossilizados é na parte superior da 
Formação. Sendo que a Formação Missão Velha pode ser segmentada em duas frações, 
conforme critérios estratigráficos. É oportuno frisar, que o objetivo deste estudo foi atingido, 
conforme verificado na discussão proposta. No entanto, este trabalho apresentou como 
principais limitações a restrição quanto ao acesso a artigos publicados de forma não pública. 
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RESUMO: Este trabalho explora eventos 
cruciais desde a Revolução Industrial, que 
desempenharam um papel fundamental na 
redução da energia cinética dos ventos e 
sua correlação com o aquecimento global. 
Enquanto os fatores que contribuem para o 
aumento das temperaturas são amplamente 
discutidos, a atenção aos elementos que 
colaboram diretamente para o resfriamento 
do planeta é notavelmente escassa. A 
circulação atmosférica, essencial nesse 
contexto, está sendo cada vez mais 
impactada. 

Na atmosfera, a coesão do ar, atuando como 
um corpo unificado devido à ausência de 
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vácuo, destaca a importância de obstáculos 
na circulação livre do ar, influenciando 
significativamente o movimento das massas 
de ar. 

O transporte de calor da superfície para 
camadas superiores da troposfera pelo 
ar úmido é crucial para a redução da 
temperatura terrestre. O resfriamento 
primário ocorre nos ventos alísios, onde o 
ar úmido ascendente, ao atingir grandes 
altitudes, irradia calor para o espaço, 
contribuindo para a redução da temperatura 
em áreas tropicais. 

A presença vital dos ventos na circulação 
atmosférica é comprometida pelo 
aumento da rugosidade da Terra devido à 
urbanização, resultando em turbulência que 
reduz a energia cinética dos ventos. 

A Revolução Industrial e o crescimento 
populacional acelerado desde 1800 
impulsionaram a concentração urbana, 
resultando em uma expansão vertical das 
cidades. Atualmente, 5% da superfície do 
planeta é ocupada por centros urbanos, e se 
prevê que 61% da população mundial viverá 
em áreas urbanas até 2025, contribuindo 
para o aumento desta estatística. 

A escassez de medidas que contribuam 
para a redução da interferência dos edifícios 
nos padrões de circulação do vento nas 
grandes cidades é alarmante, fazendo-se 
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necessária a condução de estudos e definição de diretrizes para orientar construções de 
maneira sustentável, promovendo abordagens mais equilibradas para o desenvolvimento 
urbano. 

A expansão urbana pode ser positivamente explorada com um planejamento eficiente na 
construção de novos edifícios, contribuindo para mitigar a turbulência dos ventos causada 
pela rugosidade superficial, resultando em uma maior transferência de calor da superfície 
terrestre para o espaço. 

PALAVRAS-CHAVE: Circulação atmosférica; resfriamento do planeta; expansão vertical das 
cidades; rugosidade da superfície terrestre; desenvolvimento urbano 


CONTEXTUALIZAÇÃO E DISCUSSÃO 


O ar exibe uma baixa condutividade térmica, sendo de 0,0026 w/mk, conforme 
indicado por Nick Connor, 2020. Entretanto, quando tal elemento encontra-se em 
movimento, torna-se um agente altamente influente na rápida modificação da temperatura 
de toda uma região ou até mesmo de um país como todo. Fenômenos meteorológicos, 
como as frentes frias, exemplificam este impacto, resultando em diminuições substanciais 
de até 15ºC em questão de horas sobre uma extensa área geográfica. 

Os ventos Alísios, provenientes da região dos trópicos (AYOADE 2003) vão sendo 
aquecidos na sua trajetória em direção à região do equador meteorológico, lá recebem 
ainda mais calor e umidade, pela baixa pressão e alta temperatura seguem em direção às 
camadas mais altas da troposfera. 

O vento úmido exerce papel essencial na redução da temperatura terrestre, 
decorrente de sua capacidade de transportar o calor da superfície do planeta até às 
camadas superiores da troposfera. Neste processo, o calor é então irradiado para o espaço, 
contribuindo significativamente para o mecanismo de resfriamento global. 

O gradiente vertical térmico nesse processo é de aproximadamente -6,5ºC a cada 
1000 metros de altitude, podendo atingir temperaturas abaixo de -60ºC, de acordo com a 
U.S. Navy Weather Research Facility. Esse ar já seco e gelado, ao iniciar o movimento de 
descida na direção dos trópicos (célula de Hadley), contribui para a redução da temperatura 
nessas áreas, ilustrado na Figura 1. 


Explorando o universo das ciências exatas: Teoria e aplicações Capítulo 5 


50 


Circulação 
Geral da Atmosfera Célula Polar 


Região de baixa pressão 


Região de alta pressão 
Célula 
de Ferrel 


Célula 
de Hadley 


Cólula 
de Hadley 


Célula 
de Ferrel 


Célula Polar 


Figura 1. 


Na atmosfera, tendo-se em vista a ausência de vácuo, a redução da energia cinética 
dos ventos interfere na circulação do ar em sua totalidade. Nos séculos passados, como no 
século XIll o planeta Terra se distinguia pelo predomínio de uma extensa cobertura vegetal, 
em marcante contraste com a presente realidade constatada nas áreas urbanas, conforme 
apresentado na Tabela 1 abaixo. 


Descrição do terreno Zo (mm) 
Terreno muito liso, gelo, lama 0,01 
Mar aberto em condição de calmaria 0,20 
Mar agitado 0,50 
Neve 3,00 
Gramado 8,00 
Pasto rugoso 10,00 
Campo de pouso 30,00 
Plantações 50,00 
Poucas árvores 100,00 
Muitas árvores, poucas construções 250,00 
Florestas altas 500,00 
Subúrbio 1.500,00 
Centros urbanos 3.000,00 


Tabela 1. Comprimento de rugosidade 


Fonte: Adaptado de Manwell (2002). 
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Na natureza intocada, observa-se uma transição gradual entre os biomas, 
caracterizada por um crescimento suave na altura das paisagens. Esse fenômeno contrasta 
com o ambiente urbano, onde edifícios altos frequentemente se erguem imediatamente 
próximos a estruturas mais baixas, conforme exemplificado na Figura 2 com fotos 
capturadas pelos autores. 


Cerrado próximo a Uberlândia-MG Bairro Morada da Colina, Uberlândia-MG 


Um dos principais catalisadores para o aumento da rugosidade do planeta foi a 
Revolução Industrial. Esse período crucial, marcado pela intensificação da produção e 
demanda por serviços, desencadeou um significativo êxodo rural para as áreas urbanas. 

Ainda, somado a este movimento migratório, têm-se o expressivo crescimento 
populacional que, a partir de 1800, o planeta atingiu o marco de 1 bilhão de pessoas, e, 
posteriormente, aumentando exponencialmente, em 1922 atinge-se 2 bilhões de pessoas, 
4 bilhões em 1975 e, atualmente, em 2028, atinge-se a marca extraordinária de 8 bilhões 
de habitantes no mundo. 

O aumento vertiginoso da população global, aliado à tendência de concentração 
urbana, desencadearam uma demanda exponencial por edifícios, resultando em um 
cenário de expansão vertical das cidades. 

Conforme os dados apresentados pelo renomado arquiteto e urbanista Jorge 
Wilheim, estima-se, para 2025, que 61% da população mundial habitará em áreas urbanas. 
Comparando-se esta projeção com o ano de 1975, este quantitativo era de apenas 37%, 
evidenciando-se, portanto, a marcante transição demográfica. (Boletim do Instituto de 
Estudos Avançados da Universidade de São Paulo). 

No que concerne à ocupação de áreas por centros urbanos, atualmente atinge-se 
o número de 5% da superfície do planeta. Embora o nível de ocupação possa parecer 
baixo, o distanciamento do crescimento dos centros urbanos com a realização de medidas 
adequadas para minimizar o impacto destes grandes aglomerados na circulação dos ventos 
é preocupante. Um exemplo disto é que nas leis de uso e ocupação do solo das grandes 
cidades as quais pesquisamos, não encontramos uma menção sequer a algo que leve a 
uma redução da rugosidade da região urbana com relação a circulação dos ventos, os 
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critérios para aprovação de construções altas são em sua maioria com base na capacidade 
de tráfegos das vias, os afastamentos e a altura com relação a rota de aviões. 

Reflete-se, portanto, aescassez de estudos e diretrizes que busquem avaliar e orientar 
a construção de forma a reduzir a interferência nos padrões de vento, proporcionando uma 
abordagem mais sustentável e equilibrada para o desenvolvimento urbano. 

Constando-se irregularidade na ocupação da superfície do planeta, influenciada 
por diversos fatores, como o nível de desenvolvimento, além da densidade populacional 
de cada país ou região, a interferência na circulação dos ventos apresenta considerável 
variação de uma região para outra. Sendo assim, estas áreas de maior desenvolvimento 
urbano, em decorrência do aumento da rugosidade da superfície, tendem a sentir mais 
pronunciadamente os efeitos da redução do volume de ar circulante. Ou seja, o aumento na 
rugosidade implica na menor transferência de calor para o espaço, limitando o fenômeno 
de resfriamento da temperatura terrestre desta região específica. 

As figuras 3 e 4 abaixo, imagens noturnas capturadas pela Nasa, corroboram na 
compreensão não só da irregularidade da ocupação das áreas metropolitanas, como 
também do aumento da rugosidade, criando-se obstáculos na circulação dos ventos. 


Figura 3. Estados Unidos 


Imagem noturna da superfície do planeta capturadas pela ISS divulgada pela Nasa. 
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Figura 4. Continente Europeu 


Imagem noturna da superfície do planeta capturadas pela ISS divulgada pela Nasa. 


Por fim, conforme estudos apresentados durante o evento Planet Under Pressure, 
organizado pelas Nações Unidas, estes indicam que nas próximas duas décadas, as 
áreas urbanas devem expandir-se em 1,5 milhão de quilômetros quadrados. Contudo, 
esta tendência pode ser explorada positivamente por meio do planejamento meticuloso 
da construção de novos edifícios, que, se orientado e executado de forma eficiente, 
tal processo pode contribuir na mitigação da turbulência dos ventos, ocasionada pela 
rugosidade superficial, e, consequentemente, na maior transferência de calor da superfície 
terrestre para o espaço. 


CONCLUSÃO 


O presente trabalho destaca eventos cruciais desde a Revolução Industrial que 
desempenharam um papel fundamental na redução da energia cinética dos ventos, 
correlacionando-se com o aquecimento global. Embora os fatores de aumento de 
temperatura sejam amplamente discutidos, a atenção aos elementos de resfriamento 
planetário é notavelmente escassa, com a circulação atmosférica essencial sendo cada 
vez mais afetada. A presença vital dos ventos na circulação atmosférica é comprometida 
pela rugosidade da Terra devido à urbanização, contribuindo para a redução da energia 
cinética dos ventos. 
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A região dos Trópicos emerge como particularmente sensível à redução de 
energia, pois as massas de ar, ao retornarem da troposfera, apresentam menor volume e 
temperaturas médias mais elevadas. A crescente rugosidade do planeta desde a Revolução 
Industrial, que impulsionou o desenvolvimento urbano, é incontestável, resultando na 
perda de energia cinética dos ventos e, consequentemente, na redução de sua velocidade, 
impactando no volume de ar circulante na célula de Hadley ao longo do tempo. 

Afalta de medidas para reduzir a interferência dos edifícios nos padrões de vento nas 
grandes cidades é preocupante, sendo urgente conduzir estudos e estabelecer diretrizes 
para orientar construções de forma sustentável, buscando abordagens mais equilibradas 
no desenvolvimento urbano. 

Entretanto, há uma oportunidade positiva na expansão urbana: um planejamento 
eficiente na construção de novos edifícios pode mitigar a turbulência dos ventos causada 
pela rugosidade superficial, resultando em uma maior transferência de calor da terra para 
o espaço. Assim, a integração de abordagens sustentáveis na urbanização torna-se crucial 
para enfrentar os desafios emergentes, oferecendo perspectivas promissoras para um 


desenvolvimento urbano equilibrado e resiliente no contexto das mudanças climáticas. 
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